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Ribosome Flow Model on a Ring 

Alon Raveh, Yoram Zarai, Michael Margaliot, and Tamir Tuller 


Abstract —The asymmetric simple exclusion process (ASEP) is 
an important model from statistical physics describing particles 
that hop randomly from one site to the next along an ordered 
lattice of sites, but only if the next site is empty. ASEP has 
been used to model and analyze numerous multiagent systems 
with local interactions including the flow of ribosomes along the 
mRNA strand. 

In ASEP with periodic boundary conditions a particle that 
hops from the last site returns to the first one. The mean held 
approximation of this model is referred to as the ribosome flow 
model on a ring (REMR). The REMR may be used to model both 
synthetic and endogenous gene expression regimes. 

We analyze the REMR using the theory of monotone dynami¬ 
cal systems. We show that it admits a continnnm of equilibrium 
points and that every trajectory converges to an equilibrium 
point. Eurthermore, we show that it entrains to periodic transi¬ 
tion rates between the sites. We describe the implications of the 
analysis results to understanding and engineering cyclic mRNA 
translation in-vitro and in-vivo. 

Index Terms —Monotone dynamical systems, first integral, 
asymptotic stability, ribosome flow model, entrainment, asymmet¬ 
ric simple exclusion process, mean Held approximation, mRNA 
translation, cyclic mRNA. 

I. Introduction 

The asymmetric simple exclusion process (ASEP) is an 
important model in statistical mechanics 1^ . ASEP describes 
particles that hop along an ordered lattice of sites. The dynam¬ 
ics is stochastic: at each time step the particles are scanned, 
and every particle hops to the next site with some probability 
if the next site is empty. This simple exclusion principle allows 
modeling of the interaction between the particles. Note that 
in particular this prohibits overtaking between particles. 

The term “asymmetric” refers to the fact that there is 
a preferred direction of movement. When the movement is 
unidirectional, some authors use the term totally asymmetric 
simple exclusion process (TASEP). ASEP was first proposed 
in 1968 ||25]| as a model for the movement of ribosomes along 
the mRNA strand during gene translation. In this context, the 
lattice models the mRNA strand, and the particles are the 
ribosomes. Simple exclusion corresponds to the fact that a 
ribosome cannot move forward if there is another ribosome 
right in front of it. ASEP has become a paradigmatic model 
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for non-equilibrium statistical mechanics il, 10 . It is used 
as the standard model for gene translation ll52ll . and has also 
been applied to model numerous multiagent systems with 
local interactions including traffic flow, kinesin traffic, the 
movement of ants along a trail, pedestrian dynamics and ad- 
hoc communication networks EH, ES- 

In ASEP with open boundary conditions, the lattice bound¬ 
aries are open and the first and last sites are connected to two 
external particle reservoirs that drive the asymmetric flow of 
the particles along the lattice. In ASEP with periodic boundary 
conditions the lattice is closed, so that a particle that hops from 
the last site returns to the first one. In particular, the number 
of particles on the lattice is conserved. 

Recently, the mean held approximation of ASEP with open 
boundary conditions, called the ribosome flow model (REM), 
has been analyzed using tools from systems and control the¬ 
ory ED, ESI, ED, ED, ED, E3- in this paper, we consider 
the mean field approximation of ASEP with periodic boundary 
conditions. This is a set of n deterministic nonlinear first-order 
ordinary differential equations, where n is the number of sites, 
and each state-variable describes the occupancy level in one of 
the sites. We refer to this system as the ribosome flow model 
on a ring (REMR). To the best of our knowledge, this is the 
first study of the REMR using tools from systems and control 
theory. 

In the physics literature, many properties have been proven 
for the ASEP with periodic boundary conditions and homo¬ 
geneous transition rates (see, e.g. CD). Another case that is 
amenable to analysis is where the transition rates vary but 
depending on the particles rather than on the sites (see e.g. El 
and the references therein). However, an analytical understand¬ 
ing of ASEP with site-dependent inhomogeneous transition 
rates is still pending. In contrast, most of the results in this 
paper hold for the general case of an inhomogeneous REMR. 

The REMR may model, for example, the translation of 
a circular mRNA or DNA molecule. Circular RNA forms 
have been described in all domains of life (see for example, 
ifTOl . IH . H, El, ifirl . ID, 0). Speciflcally, it was shown 
that in prokaryotes it is possible to regulate translation from 
circular DNA IH, 0- In the case of eukaryotes, the canonical 
scanning model requires free ends of the mRNA ||2T1 : however, 
it is well-known that in these organisms mRNA is often 
(temporarily) circularized by translation initiation factors EOl . 

We show that the REMR admits a continuum of equilibrium 
points, and that every trajectory converges to an equilibrium 
point. Eurthermore, if the transition rates between the sites 
vary in a periodic manner, with a common period T, then every 
trajectory converges to a periodic solution with period T. In 
other words, the REMR entrains (or phase-locks) to the peri¬ 
odic excitation. In the particular case where all the transition 
rates are equal all the state variables converge to the same 
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value, namely, the average of all the initial values. We discuss 
the implications of these results to mRNA translation. 

The remainder of this paper is organized as follows. Sec¬ 
tion |II] reviews the RPMR. Section Hn] details the main results. 
To streamline the presentation, the proofs are placed in Ap¬ 
pendix A. The final section summarizes and describes several 
possible directions for further research. 

We use standard notation. For an integer i, in € K" is the n- 
dimensional column vector with all entries equal to i. For a 
matrix M, M' denotes the transpose of M. We use | • |i : 
R" —>■ R+ to denote the Li vector norm, that is, \z\i = \zi\ + 
■ • • + \zn\- For a set K, int(Ar) is the interior of K. 

II. The model 

The ribosome flow model on a ring (RFMR) is given by 

±1 = XnXnil - Xi) - Aia:i(l - X 2 ), 

±2 = AiXi(l - X 2 ) - A 2 X 2(1 - X 3 ), 

±3 = A 2 X 2(1 - X 3 ) - A 3 X 3(1 - X 4 ), 

Xn —1 — Xn— 2 Xn — 2 (^ Xn—l) Xn—lXn—l{^ 

in = A„_ia;„_i(l - Xn) - A„x„(l - xi). (I) 

Here Xi{t) S [0,1] is the normalized occupancy level at site i 
at time t, so that Xi{t) = 0 [xi{t) = 1] means that site i is com¬ 
pletely empty [full] at time t. The transition rates Ai,..., A„ 
are all strictly positive numbers. To explain this model, con¬ 
sider the equation ±2 = AiXi(l — X 2 ) — A2X2(1 — X 3 ). The 
term ri 2 := Aia:i(l —^ 2 ) represents the flow of particles from 
site 1 to site 2. This is proportional to the occupancy xi at 
site 1 and also to 1 — X 2 , i-e. the flow decreases as site 2 
becomes fuller. This is a relaxed version of simple exclusion. 
The term r 23 := X 2 X 2 {^ — X 3 ) represents the flow of particles 
from site 2 to site 3. The other equations are similar, with the 
term r„i := A„x„(l — xi) appearing both in the equations 
for xi and for due to the circular structure of the model 
(see Fig. [Til- 

The RFMR encapsulates simple exclusion, unidirectional 
movement along the ring, and the periodic boundary condition 
of ASER This is not surprising, as the RFMR is the mean held 
approximation of ASEP with periodic boundary conditions 
(see, e.g., E] p. R345] and gH p. 1919]). 

Note that we can write (ID succinctly as 

Xi = Ai_iXi_i(l - Xi) - Xixfll - Xi+i), i = 1,. . .,n, 

where here and below every index is interpreted modulo n. 
Note also that 0„ [1„] is an equilibrium point of ([D- Indeed, 
when all the sites are completely free [completely full] there 
is no movement of particles between the sites. 

Let 

C" :={yeR" [0,1], i = l,...,n}, 

i.e., the closed unit cube in R". Since the state-variables 
represent normalized occupancy levels, we always consider 
initial conditions x(0) € C". It is straightforward to verify 
that C" is an invariant set of (ID- i-e. x(0) € C" implies 
that x{t) G C" for all t > 0. 


mRNA or DNA 



Transition rate 

- ►Al A2 AZ 



Fig. 1. A. Illustration of a cyclic mRNA or DNA molecule. With this 
topology the ribosomes terminating translation may re-initiate ED, El. d 
translation with high probability. B. The RFMR as a model for translation 
with cyclic mRNA. In this context, Xi(t) is the normalized ribosome density 
at site i at time t. The transition rates Xi depend on factors such as tRNA 
abundance. 


Note that ([D implies that 

n 

'Y2 Xi{t) = 0 , for all f > 0, 

i=0 

so the total occupancy H(x) := I'nX is conserved: 

H{x{t)) = H{x{Q)), forall<>0. (2) 

The dynamics thus redistributes the particles between the sites, 
but without changing the total occupancy level. In the context 
of translation, this means that the total number of ribosomes 
on the mRNA is conserved. 

Eq. (ID means that we can reduce the n-dimensional RFMR 
to an (n — 1)-dimensional model. In particular, the RFMR 
with n = 2 can be explicitly solved. In this case we can 
obtain explicit expressions for important quantities, e.g., the 
rate of convergence to equilibrium. This solution is detailed 
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in Appendix B. 

If we change the first [last] equation in O to Xi = 
Aq (1 AiXi(l X 2 ) \_Xn — ^n — lXn — l (1 ^n) AyjXyj] 

we obtain the RFM. This may seem like a minor change, 
but in fact the dynamical properties of the RFM and the 
RFMR are very different. For example, in the RFM there is no 
first integral. Also, in the RFM there is a single equilibrium 
point in C", whereas as we shall see below the RFMR has a 
continuum of equilibrium points in C". 

The next section describes several theoretical results on 
the RFMR. Since the case n = 2 is solved in Appendix B, we 
assume from here on that n > 3. Applications of the analysis 
to gene translation are discussed in Section HVl 

III. Main results 
A. Strong Monotonicity 

A cone K C K” defines a partial ordering in K” as follows. 
For two vectors a,b G M", we write a < b if {b — a) G 

K', a < 6 if a < b and a ^ b', and a <C b if {b — a) G 

int(iT). The system y = f{y) is called monotone if a < b 
implies that y{t, a) < y{t, b) for all f > 0. In other words, the 
flow preserves the partial ordering ll42l . It is called strongly 
monotone if a < 6 implies that y{t, a) <C y{t, b) for all t > 0. 

From here on we consider the particular case where the 
cone is K = R". Then a < 6 if < 6^ for all i, and a b 

if Qi < bi for all i. A system that is monotone with respect to 

this partial ordering is called cooperative. 


Proposition 1 Let x{t, a) denote the solution of the RFMR at 
time t for the initial condition x(0) = a. For any a,b G C" 
with a < b we have 


■+-r 

VI 

■+-r 

for all t > 0. 

(3) 

Furthermore, if a < b then 



x{t, a) <C x{t, b), 

for all t > 0. 

(4) 


In the context of translation, this means the following. Con¬ 
sider two possible initial ribosome densities on the same cyclic 
mRNA strand, a and b with < bi for all i, that is, b 

corresponds to a higher ribosome density at each site. Then 
the trajectories x{t, a) and x{t, b) emanating from these initial 
conditions continue to satisfy the same relationship between 
the densities for all time t >0. 

B. Stability 

Denote the s level set of FI by 

:= {y € C" : I'^y = s}. 

The next result shows that every level set contains a unique 
equilibrium point, and that any trajectory of the RFMR em¬ 
anating from any point in Ls converges to this equilibrium 
point. In the context of translation on a cyclic mRNA strand, 
this means that a perturbations in the distribution of ribosomes 
along the strand (that does not change the total number of 
ribosomes) will not change the asymptotic behavior of the 
dynamics. It will still converge to the same unique steady state 



Fig. 2. Trajectories of o with n = 3 for three different initial conditions 
in L 2 '. [1 1 0]', [1 0 1]', and [0 1 l]f The equilibrium point is marked 
with a circle. 


ribosome distribution and therefore lead to the same steady- 
state translation rate. 


Theorem 1 Pick s G [0, n]. Then Lg contains a unique 
equilibrium point of the RFMR, and for any a G Lg, 

lim x{t, a) = Cl,. 

i—>oo 

Furthermore, for any 0 < s < p < n, we have 

eL 3 <eip. (5) 

Thm. [T] implies that the RFMR has a continuum of linearly 
ordered equilibrium points, namely, {eL„ ■ s G [0, n]}, 
and also that every solution of the RFMR converges to an 
equilibrium point. 


Example 1 Consider the RFMR with n = 3, Ai = 2, A 2 = 3, 
and A3 = 1. Fig. |2] depicts trajectories of this RFMR for three 
initial conditions in L 2 : [1 1 0]', [1 0 1]', and [0 1 1]'. It 
may be observed that all the trajectories converge to the same 
equilibrium point 6^2 ~ [0.5380 0.6528 0.8091] ^ Fig. |3] 
depicts all the equilibrium points of this RFMR. Since A 2 > 
Ai and A2 > A3, the transition rate into site 3 is relatively 
large. As may be observed from the figure this leads to 63 > ei 
and 63 > 62 for every equilibrium point e. □ 

Fix an arbitrary s G [0, n]. To simplify the notation, we just 
write e instead of from here on. Then 

l^e = s, 


( 6 ) 
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Fig. 3. All the equilibrium points of the RFMR with n = 3, Ai =2, 
A 2 = 3, and A 3 = 1. 


and since for x = e the left-hand side of all the equations 
in o is zero, 

Arie„(l - ei) = Aiei(l - 62 ) 

= A2e2(l — 63) 

— — 671). ( 7 ) 

Thus, the flow along the chain always converges to a steady- 
state value 

R := Xiei{l — Ci+i), for all i. 

Using (|7]i and the equation l^e = s it is possible to derive 
a polynomial equation for ei. For example, when n = 2 

(A 2 — Ai)e^ -f (Ai(s — 1) — A 2 (s -f l))ei -I- A 2 S = 0, 

whereas for n = 3 we get 

A 1 A 3 (Ai -f A 2 -f A 3 ) ef 

+ (A 2 A 3 — Xi (A 2 + A 3 S) — A 3 A 1 (A 2 ( 2 s — 1) -f A 3 (s + 1))) e 
+ (A 3 A 1 (A 2 (s^ — 3 ) -f A3(2 s — 1)) 

+ Aj (A2 (s — 2) + Xs^s — 1)) — A2A3(s + ‘2))ei 
+ A 3 (A 2 A 3 ( 2 s -f 1 ) -f Ai (A 2(1 — (s — 2 )s) — A 3 (s — 1 ))) ei 
— A2A3S = 0. 

These equations can be solved numerically, but their analysis 
seems non trivial. 

C. Differential Analysis 

Differential analysis is a powerful tool for analyzing non¬ 
linear dynamical systems (see, e.g., m, ESI, ID). The basic 
idea is to study the difference between trajectories that em¬ 
anate from different initial conditions. The next result shows 
that the RFMR is non-expanding with respect to the Li norm. 


Proposition 2 For any a, 6 € C”, 

\x{t,a) — x(t,b)\i < \a — b\i, forallt>0. ( 8 ) 

In other words, the Li distance between trajectories can never 
increase. From the biophysical point of view, this means that 
the Li difference between two profiles of ribosome densities, 
related to two different initial conditions, in the same cyclic 
mRNA/DNA is a non-increasing function of time. 

Example 2 Pick a,b £ C” such that 6 < a. By monotonic¬ 
ity, x{t,b) < x{t,a) for all f > 0 , so d{t) := |x(f, a) — 
x{t, &)|i = ln(x{t, a) — x{t, b)). Thus, 

d{t) = l'^x{t,a) - lnx{t,b) 

= 0 - 0 , 

so clearly in this case ® holds with an equality. □ 

Pick a £ C”, and let s := Substituting b = in ® 
yields 

a) - eiji < |a-Ciji, for alH > 0. (9) 

This means that the convergence to is monotone in 
the sense that the Li distance to can never increase. 
Combining (|9]l with Theorem [T] implies that every equilibrium 
point of the RFMR is semistable M- 

D. Entrainment 

Suppose now that the transition rates along the cyclic 
mRNA (or DNA) molecule are periodically time-varying 
functions of time with a common (minimal) period T > 0. 
This may correspond for example to periodically varying 
abundances of tRNA due to the cell-division cycle that is a 
periodic program for cell replication. A natural question is will 
the mRNA densities along the mRNA strand (and thus the 
translation rate) converge to a periodically-varying pattern? 

We can study this question using the RFMR as follows. We 
say that a function / is T-periodic if /(f-|-T) = f{t) for all t. 
Assume that the A^s are time-varying functions satisfying: 

• there exist 0 < i5i < (52 such that Xi{t) £ [ 5 i,( 52 ] for 
all f > 0 and alH G {1 ,..., n}. 

• there exists a (minimal) T > 0 such that all the A^s are T- 
periodic. 

We refer to the model in this case as the periodic ribosome 
flow model on a ring (PRFMR). 

Theorem 2 Consider the PRFMR. Fix an arbitrary s £ [0, n]. 
There exists a unique function fs '■ K-i- C”". that is T- 
periodic, and 

lim a) — 0 s(f)| = 0 , foralla£Ls. 

t—^OO 

In other words, every level set Lg of H contains a unique 
periodic solution, and every solution of the PRFMR emanating 
from Lg converges to this solution. Thus, the PRFMR entrains 
(or phase locks) to the periodic excitation in the A^s. 

Note that since a constant function is a periodic function 
for any T, Thm.|2] implies entrainment to a periodic trajectory 
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Fig. 4. Solution of the PRFMR in Example[^ solid line-xi (f, a); dash-dotted 
line-X 2 (t, a); dotted line-X 3 (t, a). 


in the particular case where one of the A^s oscillates, and all 
the other are constant. Note also that Thm. [T] follows from 
Thm.|2] 


Example 3 Consider the RFMR with n = 3, Xi{t) = 3, 
A 2 (f) = 3 + 2sin(f + 1/2), and Xsit) = 4 — 2cos(2f). 
Note that all the A^s are periodic with a minimal com¬ 
mon period T = 27r. Fig. |4] shows the solution x{t, a) 
for a = [0.50 0.01 0.90]/ It may be seen that every Xi(t) 
converges to a periodic function with period 27r. □ 

All the results above hold for general rates A^. In the 
particular case where all the rates are equal, it is possible to 
provide stronger results. 


E. The homogeneous case 

It has been shown experimentally that in some cases the 
ribosomal elongation speeds along the mRNA sequence are 
approximately equal El. To model this case, assume that 

Ai — * * * — Xji .— Ac, 

i.e., all the transition rates are equal, with Ac denoting their 
common value. In this case O becomes: 

Xi = XcXn(l - Xi) - XcXi{l - X2), 

±2 = XcXi(l - X2) - AcX 2(1 - X3), 

Xn — AcXyj_i(l Xn^ XcXyi(^\ ^l')’ (10) 


We refer to this as the homogeneous ribosome flow model on 
a ring (HRFMR). Also, © becomes 

e„(l - ei) = ei(l - 62 ) 

= 62(1 - 63) 

— f-n— 1(1 f-n), (11) 

and it is straightforward to verify that e = cl„, c £ K, 
satishes E). 

Dehne the averaging operator Ave(-) : R" —>■ R 
by Ave( 2 ) := il/ 2 . 

Corollary 1 For any a € C” the solution of the HRFMR 
satisfies 

lim x{t,a) = Ave(a)l„. (12) 

From the biophysical point of view, this means that equal 
transition rates along the circular mRNA lead to convergence 
to a uniform distribution of the ribosome densities. 

Note that (fT2l i implies that the steady-state flow is i? = 
Ac Ave(a)(l—Ave(a)). Thus, i? is maximized when Ave(a) = 
1/2 and the maximal value is R* = Ac/ 4 . 

Remark 1 It is possible also to give a simple and self- 
contained proof of Corollary [T| using standard tools from the 
literature on consensus networks J30l. Indeed, pick r > 0 and 
let i be an index such that xflr) > Xj{T) for all j ^ i. Then 

xflr) = Xi-i{T){l - xflr)) - Xi{T){l - Xi+i{T)) 

< Xi{T){l - Xiir)) - Xi{T){l - xflr)) 

= 0 . 

Furthermore, if Xi{T) > Xj{T) for all j i then Xi{T) < 0. A 
similar argument shows that if ^^(t) < Xj{T) {xflr) < Xj{T)] 
for all j i then Xi{T) > 0 {xflr) > 0]. Thus, the max¬ 
imal density never increases, and the minimal density never 
decreases. Dehne V{-) : M" —>■ M+ by V{y) := max^ — 
minimi. Then V{x{t)) strictly decreases along trajectories of 
the HRFMR unless x{f) = cl„ for some c £ R, and a standard 
argument (see, e.g., 1231 ) implies that the system converges to 
consensus. Combining this with © completes the proof of 
Corollary [T] 

We note in passing that Corollary [T]implies that the HRFMR 
may be interpreted as a nonlinear average consensus network. 
Indeed, every state-variable replaces information with its two 
nearest neighbors on the ring only, yet the dynamics guarantee 
that every state-variable converges to Ave(a). Consensus net¬ 
works are recently attracting considerable interest l30l . 1^ . 
1 ^ . and have many applications in distributed and multi-agent 
systems. 

The physical nature of the underlying model provides a 
simple explanation for convergence to average consensus. 
Indeed, the HRFMR may be interpreted as a system of n 
water tanks connected in a circular topology through identical 
pipes. The How in this system is driven by the imbalance in the 
water levels, and the state always converges to a homogeneous 
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distribution of water in the tanks. Since the system is closed, 
this corresponds to average consensus. 

We next analyze the linearized model of the HRFMR near 
an equilibrium point to obtain information on the convergence 
rate and the amplitude of the oscillations. 

1) Local analysis: Every trajectory of the HRFMR con¬ 
verges to cl„, where c depends on the initial condition. 
Let y := x — cl„. Then a calculation shows that the linearized 
dynamics of y is given by 

y = Qy, ( 13 ) 


where 


-1 c 0 0 

1 -c -1 c 0 

0 1 -c -1 c 


0 1 -c 

0 0 

0 0 


(14) 



_c 0 0 0 ... l-c-l_ 

By known-results on circulant matrices (see, e.g., El), the 
eigenvalues of Q are 


Fig. 5. log(|x(t) — ( 1 / 4 ) 14 !) iri the HRFMR with n = 4 and 3:(0) = 
[1 0 0 0]' as a function of t (solid line). Also shown is the function —t + 
log(^3/4) that is obtained from the local analysis (dashed line). 


7 ^ = -1-l-cfc^ ^-f (1 — c)fy^^ f=l...,n, (15) 

where w := exp(27rV—1/n), and the corresponding eigen¬ 
vectors are 


the graph of — f -f \og{^J2>/4). It may be seen that the real 
convergence rate is slightly faster than the estimate ( fTSl l. □ 


/:=[! ... £=l...,n. 

(16) 

In particular, 71 = 0, with corresponding eigenvector = 1„. 
This is a consequence of the continuum of equilibria in the 
HRFMR. 

Note that 

Re( 7 ^) = —1 -j- cos(27r(^ — l)(n — l)/n) 

-I- c(cos(27r(f — l)/n) — cos(27r(£ — l)(n — l)/n)) 
= —1-I-cos(27r(f — l)/n), (17) 

and this implies that 

Re( 7 f) < Re( 72 ) = cos(27r/n) — 1, i = 2,... ,n. 

Thus, for x(0) in the vicinity of the equilibrium 

\x(t) — cl„| < exp((cos(27r/n) — l)f)|x(0) — cl„|. (18) 

The exponential convergence rate decreases with n. For ex¬ 
ample, for n = 2, cos(27r/n) — 1 = —2, whereas for n = 10, 
cos(27r/n) — 1 Ri —0.191. In other words, as the length of 
the chain increases the convergence rate decreases. This is the 
price paid for the fact that each site “communicates” directly 
with its two neighboring sites only. 

Our simulations suggest that (fTsT i actually provides a rea¬ 
sonable approximation for the real convergence rate (i.e., not 
only in the vicinity of the equilibrium point) in the HRFMR. 
The next example demonstrates this. 

Example 4 Consider the HRFMR with n = 4. Fig. |5] de¬ 
picts log(|x(f) — (1/4)14!) for rhe initial condition x(0) = 
[100 O]^ Note that here log(|x(0) — (l/4)l4|) = 
log( 1 / 3 / 4 ). In this case, Re( 72 ) = —1, so (fTSl l becomes 
log(|x(f) — cl„|) Ki —t + log(|x(0) — cl„|). Also shown is 


Eq. (fTTl i implies that c does not affect the convergence rate 
in the linearized system. It does however affect the oscillatory 
behavior until convergence. To see this, note that the solution 
of (fTsT l is 

n 

y(t) = ^ Sz exp( 7 if)z;/ (19) 

where the satisfy 

n 

y(o) = 

i=l 

Since 71 = 0 and Re( 7 i) < 0 for all i > 1, ( fT9l l implies 
that limi_,.oo y(i) = so si = 0. The three eigenvalues 

with the largest real part are 71 , 72 , and 7 „, so for large t, 

y{t) R! S 2 exp( 72 f)v^ -|- s„ exp( 7 „f)w”. 

It follows from (flSl l that 

a := Re( 72 ) = Re( 7 „) = cos(27r/n) — 1, 

/? := Im( 72 ) = - Im( 7 „) = (2c - 1) sin(27r/n), 

so u" = if^, Sn = S 2 and this yields 

y{t) R 2|p| exp(af) cos (/?f + Zp), 

where p := S 2 V^. The oscillatory behavior thus depends on c. 
For c = 1/2 the solution has no oscillations at all, and as c 
moves away from 1/2 the oscillations become larger. The 
reason for this is that for c = 1/2 the linear equation (fOl l 
corresponds to the case where each agent weighs the con¬ 
tribution from its two neighbors equally. As c moves away 
from 1/2 the weights become different and this imbalance 
leads to oscillations until the state-variables converge to the 
correct values. 
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Fig. 6. The function d{t, c) for several values of c. 


Examples Let z(c) := [c + 0.1 c—0.1 c\'. We simu¬ 

lated the trajectories of the HRFMR with n = 3 for four 
initial conditions: Xq = z(c), with c = 0.1,0.3,0.5, and 0.9. 
Note that ^l'^z{c) = c, so limt_,.oo a:(f, z(c)) = CI 3 . Fig. | 6 ] 
depicts 

d(f,c) := {xi{t,z{c)) - c) - (xi(f,z(l/ 2 )) - 1 / 2 ) 

as a function of t. It may be seen that as |c — 1/2| increases 
the oscillations increase. This agrees with the analysis above. 


E The Case of a Single Slow Rate 

It is interesting to consider the case where all the transition 
rates are equal to Ac, except for Ai that has value Xg, 
with Xq < Ac. For TASEP with periodic boundary conditions 
this case has been studied in 03, HI. It corresponds to a 
translation regime with the elongation rates basically uniform 
(with rate Ac), yet the re-initiation, and possibly termination, 
rate (modeled by Xq) is slower. In this case, Eq. (|2]l of the 
REMR becomes 

en(l - ei) = ^ei(l - 62 ) 

= 62(1 - eg) 

= e„_i(l-e„), ( 20 ) 

so we assume from here on, without loss of generality, 
that Ac = 1 . 

When 1 ^ — 5 1 is small, i.e., the normalized total occupancy is 
close to 1/2 the steady-state e has the form depicted in Eig.|2] 
The slow rate yields an increase [decrease] in the steady-state 
particle density to the immediate left [right]. In other words, 
the slow rate induces a “traffic jam” segregating the steady- 
states into high- and low-density regions. Near the slow site the 
densities become more or less uniform with a low density e/ 



Fig. 7. Steady-state occupancy levels e^, z = 1,.. ., 40, in a RFMR with 
n = 40, s = 20.3 and A 2 = ■ ■ ■ = A^, = 1. Upper figure: Ai = 0.05 
(’+’) and Ai = 1 (’o’) (i.e., the HRFMR). Lower figure: Ai = 0.5 (’+’) and 
Ai = 1 (’o’). 


and a high density eh- Substituting this in (l20l) gives 

e/i(l - eh) ^ Xqehil - ei) 

K. eiil - ei), 

so ei ~ and eh ~ Note that this implies that e; -f 

eh ~ 1. Eor example, for Xq = 0.05 this gives ei « 0.047619, 
eh ~ 0.952381, and this agrees well with the case depicted in 

Fig. ID 

The steady-state flow is thus 

R = eh{l - eh) 

... A, 

^ (1 + A,) 2 - 

Let mi [mh] denote the number of e^s satisfying d « ei, 
so that mh := n — mi approximates the number of e^s 
satisfying ~ eh- Then the equation s ~ miei + mheh yields 

n — s(l -I- Xq) 

T^a; ■ 

Eor example, for the case n = 40, s = 20.3, Xq = 0.05 
this gives mi « 19.6684, and this agrees well with the case 
depicted in Eig. [T] 

Erom the biophysical point of view, these results suggest 
that a slower re-initiation step (that may interact and delay 
the termination step) will lead to a ribosomal “traffic jam” 
before the STOP codon. 

IV. Discussion 

The ribosome flow model on a ring (REMR) is the mean 
field approximation of ASEP with periodic boundary condi¬ 
tions. We analyzed the REMR using tools from monotone 
dynamical systems theory. Our results show that the REMR 
has several nice properties. It is an irreducible cooperative 
dynamical system admitting a continuum of linearly ordered 
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equilibrium points, and every trajectory converges to an equi¬ 
librium point. The RFMR is on the “verge of contraction”, 
and it entrains to periodic transition rates. 

Topics for further research include the following. ASEP 
with periodic boundary conditions has been studied exten¬ 
sively in the physics literature and many explicit results are 
known. For example, the time scale until the system relaxes to 
the (stochastic) steady state is known 0 . A natural research 
direction is based on extending such results to the RFMR. 

For the RFM, that is, the mean-field approximation of ASEP 
with open boundary conditions, it has been shown that the 
steady-state translation rate R satisfies the equation 

0 = /(i?), 

where / is a continued fraction ll28l . Using the well-known 
relationship between continued fractions and tridiagonal ma¬ 
trices (see, e.g., Il49ll ) yields that R~^l'^ is the Perron root of a 
certain non-negative symmetric tridiagonal matrix with entries 
that depend on the A^s llTSl . This has many applications. For 
example it implies that i? = i?(Ao,..., A^) in the RFM is a 
strictly concave function on llTSl . It also implies that 

sensitivity analysis in the RFM is an eigenvalue sensitivity 
problem 1^ . An interesting research question is whether R 
in the RFMR can also be described using such equations. 

Another interesting topic for further research is studying a 
network of several connected RFMs. The output of each RFM 
is divided between the inputs of all the RFMs. This models 
competition for the ribosomes in the cell. Assuming that 
the system is closed then leads to network of intercon¬ 
nected RFMRs. 

The RFMR may be used in the future for analyzing 
novel synthetic circular DNA or mRNA molecules for protein 
translation in-vitro or in-vivo. Such devices have potential 
advantages with respect to linear mRNA, since they do not 
include free ends and may thus be more stable. For example, 
in E. coli RNA degradation often begins with conversion of 
the 5'-terminal triphosphate to a monophosphate, creating a 
better substrate for internal cleavage by RNase E iJTl . In 
eukaryotes, such as Saccharomyces cerevisiae, intrinsic mRNA 
decay initiate with deadenylation that causes the shortening of 
the poly(A) tail at the 3' end of the mRNA, followed by the 
removal of the cap at the 5' end by the decapping enzyme, 
which leads to a rapid 5' —3' degradation of the mRNA by 
an exoribonuclease ioi. In eukaryotes, the canonical scanning 
model requires free ends of the mRNA ET\ . However, it may 
be possible to design circular DNA or mRNA molecules by 
using Internal Ribosome Entry Sites (IRESes) HTII . Such an 
experimental system (or a similar system) can be used in the 
future for evaluation the theoretical results reported in this 
study. 

Appendix A; Proofs 

Proof of Prop. Write the RFMR ([T]i as i = f{x). The 
Jacobian matrix J{x) := ^(x) is given in ( 1211 1 . This matrix 
has nonnegative off-diagonal entries for all x G C". Thus, 
the RFMR is a cooperative system ll42l . and this implies ([^l. 
Furthermore, it is straightforward to verify that J(x) is an 


irreducible matrix for all x € 101(17"), and this implies (|4]i 
(see, e.g., ||42l Ch. 4]). ■ 

Proof of Thm. Q] Since the RFMR is a cooperative irre¬ 
ducible system with H{x) = V^x as a first integral, Thm. [T] 
follows from the results in (see also ED and ES for 
some related ideas). ■ 

Proof of Prop. |2l Recall that the matrix measure pi{-) : 
]R"x™ induced by the Li norm is 

p.i{A) = max{ci(A),. .. ,Cn(A)}, 

where Ci(A) := an + Wki\, i-e. the sum of entries 

in column i of A, with the off-diagonal entries taken with 
absolute value ll48l . For the Jacobian of the RFMR, we 
have Ci{J{x)) = 0 for all i and all x € (7", so pi{J{x)) = 0. 
Now (Is) follows from standard results in contraction theory 
(see, e.g., EH)- ■ 

Proof of Thm. |2l Write the PRFMR as i = f{t,x). 
Then f{t,y) = f(t + T,y) for all t and y. Further¬ 
more, H{x) = ll^x is a first integral of the PRFMR. Now 
Thm. |2] follows from the results in ED (see also GQl). ■ 
Proof of Corollary [J Fet s := V^a. Then Eg con¬ 
tains Ave(a)l„ and this is an equilibrium point. The proof 
now follows immediately from Thm. [D ■ 

Appendix B: Solution of the RFMR with n = 2 
Consider ([D with n = 2 , i.e. 

±1 = \2X2{1 — xi) — Aia;i(l — X 2 ), 

±2 = XlXi{l — X 2 ) — X2X2{i — Xi). ( 22 ) 

We assume that a;(0) 7 ^ O 2 and a;(0) 7 ^ I 2 , as these are 
equilibrium points of the dynamics. Fet s := a;i(0) -|- a; 2 ( 0 ). 
Substituting X 2 (t) = s — Xi{t) in (l22l i yields 

Xi = A 2 (s — a;i)(l — Xi) — AiXi(l — s -f Xi) 

= a2xf + aixi + ao, (23) 

where 

CH2 ■= A 2 — Ai, 

O'! := (Ai — A 2 )s — Ai — A 2 , 
ao •= SA 2 . 

If Ai = A 2 then (l23t is a linear differential equation and its 
solution is 

xi{t) = ^(1 — exp(—2Aif)) -I- a;i(0) exp(—2Aif), (24) 
so 

X2{t) = S — Xi{t) 

= ^(1 + exp(— 2 Aif)) — a;i( 0 ) exp(— 2 Aif) 

= |(1 - exp(-2Aif)) -I- a; 2 ( 0 ) exp(-2Aif). (25) 
In particular, 

lim x{t) = (s/ 2 )l 2 , (26) 

t—^OO 

i.e., the state-variables converge at an exponential rate to the 
average of their initial values. 


0 

0 

0 


An(l - Xx) 
0 
0 
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J{x) = 


-XnXn — Ai(l — X 2 ) 
Ai(l - X2) 

0 


\lXi 

-Aixi - A2(1 - X3) 
A 2(1 - 2 : 3 ) 

0 

0 


0 

A 2 a :2 

— A2X2 — A 3(1 — X4) 

0 

0 


— Xn.-2Xn.-2 — — Xn) Xn-lXn -1 

An—1(1 ~ x:n) —Xn—lXn —1 — An(l — Xi) 


( 21 ) 


In the context of translation on a circular mRNA/DNA this 
means that uniform translation rates are expected to yield 
a steady-state of uniform ribosomal distributions along the 
transcript, and that the convergence to this steady-state is fast. 


If Ai 7 ^ A 2 then ( |2^ is a Riccati equation (see, e.g. IfT^ l. 
whose solution is 


Xi{t) 


where 


— (y.\ — A/Acoth(-\/A(f — fo)/2) 
2 a 2 


(27) 


A := — 40 : 20:0 

2 , -1 

to := —coth 
VA 


— (s — l)^(Ai — A 2 )^ + 4 A 1 A 2 , 
/ 2x1(0)02 -f oi A 

I VK )■ 


Note that since the A^s are positive, A > 0. Also, a 
straightforward calculation shows that tg is well-defined for 
all xi(0) G [0,1]. Note that dZTl i implies that 

lim x{t) = -— T—Oi — a/A 2 a 2 S + ai + '. 

t^oo 2 o 2 


The identity 


coth f ^a/A ) — 1 =-i- (28) 

V 2 J exp(yAt)-l 

implies that for sufficiently large values of t the convergence 
is with rate exp(—v^t). Thus, the convergence rate depends 
on Ai, A 2 , and s. 



Fig. 8. Trajectories of ffl with n = 2, Ai = 2 and A 2 = 1 for three initial 
conditions. The dynamics admits a continuum of equilibrium points marked 
by -I-. 


where 


and 


\/3 + (s-l )2 
- 2 -’ 

k := tanh“^ ((xi(0) -f 1 — s/2)/z). 


Summarizing, when n = 2 every trajectory of the RFMR 
follows the straight line from x(0) to an equilibrium point e = 
e(Ai, A 2 , s). In particular, if a, 6 G satisfy l^a = I 2 & then 
the solutions emanating from a and from b converge to the 
same equilibrium point. Fig. [8] depicts the trajectories of the 
RFMR with n = 2, Ai = 2 and A 2 = 1 for three initial 
conditions. 


To study entrainment in this case, consider the RFMR 
with n = 2, Ai(f) = 3q{t)/2, and A 2 (f) = q{t)/2, where q{t) 
is a strictly positive and periodic function. Then (|2^ becomes 

X 1 = {—xf -f (s — 2)xi -I- s/2)q. (29) 

Assume that 

xl{0)<s/2. (30) 


It is straightforward to verify that in this case the solution 

of (I29I 1 is 


xi{t) = (s/2) -1-1- ztanh Ik + z 


q{s)ds , 


Jo 


Note that (l30l l implies that k is well-defined. Suppose, for 
example, that q{t) = 2 + sin(f). Then Ai(f) and A 2 (t) are 
periodic with period T = 2tt. In this case, 

xi{t) = (s/2) -1-1- 2 ;tanh (k + z(2t -f 1 — cos(f))), 

and 


X2{t) = S - Xi(f) 

= (s/2) -I- 1 — 2 ; tanh (k + z(2t -f 1 — cos(f))). 

Thus, for every a G Lg, limt_^oo a:(f, a) = where 

(f)s[t) = \[s/2) — \ + z (s/ 2 ) -I- 1 — z\ (which is of course 
periodic with period T). 
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